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Nano-pillar 


Nano-contact 


Passage  of  an  electric  current  through  a  multilayered  magnetic  nano-structure  can  lead  to  the 
excitation  of  a  persistent  magnetization  precession  in  the  thinner  ("free")  magnetic  layer  of  the 
structure.  The  frequency  of  the  excited  auto-oscillations  is  close  to  the  ferromagnetic  resonance 
frequency  and  typically  lies  in  the  microwave  frequency  range.  This  effect  can  be  used  for  the 
development  of  a  novel  class  of  fully  metallic  microwave  nano-sized  oscillators  -  spin-torque 
oscillators  (STO). 
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Equation  of  motion 
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Landau-Lifshits-Gilbert-Slonczewski  equation: 

dM 

-  =  r[HeffxM]  +  TG  +  Ts 
at 


r[Heff  X  M]  -  conservative  torque  (precession) 


damping 


precession 


Tg  = - — [Mx[MxHff]]  -  dissipative  torque  (positive  damping) 


M 

<jI 


o 


Ts  =  + - [Mx[Mx  p]]  -  spin-transfer  torque  (negative  damping) 
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STO  and  conventional  oscillator 
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Nonlinear  oscillator  model 
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oscillation 

(precession) 


positive  damping 
(Gilbert  torque) 


negative  damping 
(spin-transfer 
torque) 


dt 


+  ia>(p)c  +  r+(p)c  -  f _(p)c  =  0 


P  = 


c 


-  (dimensionless)  oscillation  power 


(/)  =  arg(c)  -  oscillation  phase 


Nonlinear  oscillator  model 


6 


Tnncusc 


> 

Oakland 

UNIVERSITY 


dc 

dt 


+  ia>{p)c  +  V+(p)c  -  r _{p)c  =  0 


Conditions  of  applicability: 


Excitation  of  only  one  mode  •  Weakly  non-conservative  system 


TOraec 


Where  is  the  magnetic  field? 


dc 

dt 


+  ico{p)c  +  TAp)c  -  T_(p)c  =  0 


All  physical  parameters  of  spin-torque  oscillators 
are  hidden  in  the  "material"  functions 
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r Ap) 
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Weakly-nonlinear  expansion: 

co(p)  =  o\  +  Np  Y+(p)  =  rG(l  +  gp)  T_(p)  =  cd(\-p) 
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Deterministic  theory  of  an  auto-oscillator 


— i-  ico(p)c  +  r+(p)c  -  r _(p)c  =  o 

dt 

Stationary  solution:  c(t )  =  exp {-i(Ogt  + 

Stationary  power  is  determined  from  the  condition  of  vanishing  total  damping: 

r+CPo)  =  r_  (p0 ) 

Stationary  frequency  is  determined  by  the  oscillation  power:  (O  =  <X>(.p$) 


Weakly-nonlinear  expansion: 


co(p)  =  co0  +  Np 

C-1 


T+(p)  =  rG(\+Qp) 


Po  = 


C  +  Q 


ox.  =  ox  +  N 
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C  +  Q 


r_(p)  =  al(l-p) 
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Stochastic  nonlinear  oscillator  model 
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Stochastic  Langevin  equation: 

^ + ia>(p)c  +  r+(p)c  -  r  (p)c  =  fn(t ) 

at 

Random  thermal  noise: 


(/,(')/;«’)) = 2  dj(i  -  o 

DJp)  =  r ApMp)  =  r,(j>)  k‘P 

Aa>(p) 


V 

Aco(p) 


-  thermal  equilibrium  power  of  oscillations: 
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Linearized  power-phase  equations 


12 


TnnzD&c 


dc 

Stochastic  Langevin  equation:  —  +  ia>{p)c  +  T+{p)c-T_(p)c  =  fn(t) 


dt 


Power-phase  ansatz:  c(t )  =  ^jp0  +  Sp(t)  e 


Sp  |«  p0 


Stochastic  power-phase  equations: 

^7  +  2reff  Sp  =  2 *Jp~Q R e[fn(t)e  ^] 

dt 

^  +  oAp„)  =  - —  Im|/,(0c-"]+  N  Sp 
dt  VA 


Linear  system  of  equations. 

Can  be  solved  in  a  general  case. 


Effective  damping:  reff  =  (G+  -G_)p0  G±=dT±(p)/dp 

Nonlinear  frequency  shift  coefficient:  N  =  dco{p)ldp 
Nonlinear  frequency  shift  creates  additional  source  of  the  phase  noise 


Linewidth  of  a  nonlinear  auto-oscillator  13 


TnnzDEC 


Lorentzian  lineshape  with  the  full  linewidth 


2A^=(i+v2)r \  k*T 


E(Po) 


=  (l  +  v2)2Aqin 


[J.-V.  Kim  etal.,  Phys.  Rev.  Lett.  100,  017207  (2008)] 


Frequency  nonlinearity  broadens  linewidth  by  the  factor 
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Linewidth  of  a  spin-torque  oscillator 
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Due  to  small  nano-scale  sizes,  spin-torque  oscillators  are  very  vulnerable  to  the  influence  of 
thermal  fluctuations.  Thermal  noise  determines  one  of  the  main,  from  the  practical  point  of 
view,  characteristics  of  STO  -  generation  linewidth. 


The  problem  of  stochastic  dynamics  of  STO  under  the  action  of  thermal  noise  was  analyzed  in 
[J.-V.  Kim,  V.  Tiberkevich,  and  A.N.  Slavin,  Phys.  Rev.  Lett.  100,  017207  (2008)],  yielding  the 
following  expression  for  the  generation  linewidth: 


2A®=r0^ 

F 


1  + 


f N ^ 

2 

{  Gj 

Here: 

r0  -  linear  Gilbert  damping  rate  (half-linewidth  of  ferromagnetic  resonance) 

kB-  Boltzmann  constant 

T- absolute  temperature 

E0  -  energy  of  auto-oscillations 

N-  nonlinear  frequency  shift  coefficient 

(rate  at  which  oscillator  frequency  changes  with  the  oscillation  power) 
G-  nonlinear  damping  coefficient 

(rate  at  which  effective  STO  damping  changes  with  the  oscillation  power) 
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The  linear  damping  rate  ro  and  nonlinear  frequency  shift  coefficient  TV for  spin-torque  oscillator 
can  be  calculated  using  standard  formulas  for  ferromegnetic  resonance.  The  energy  of  auto¬ 
oscillations  E0  and  nonlinear  damping  coefficient  G  are  given  by 


(-i 

r  (  +  Q 


G  =  ((  +  QK 


Here  M0  is  the  saturation  magnetization  of  the  "free"  magnetic  layer,  Veff  is  the  effective  volume 
of  this  layer,  involved  in  the  microwave  precession,  /is  the  modulus  of  gyromagnetic  ratio,  Q  is 
the  phenomenological  nonlinear  damping  parameter  (dimensionless  quantity  of  the  order  of 
unity),  and  <^is  the  supercriticality  of  the  bias  current: 


2eM0VQ{{ 

I  is  the  bias  current  traversing  the  magnetic  structure,  7th  =  cr/ro  is  the  threshold  value  of  this 
current,  at  which  self-sustained  auto-oscillations  starts,  £is  the  dimensionless  spin-polarization 
efficiency,  g  is  the  spectroscopic  Lande  factor,  juB  is  the  Bohr  magneton,  and  e  is  the  modulus  of 
electron  charge. 
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Nonlinear  broadening  of  STO  linewidth 
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In  contrast  with  the  majority  of  conventional  auto-oscillators,  STO  is  strongly  nonlinear 
oscillator,  in  the  sense  that  the  precession  frequency  strongly  depends  on  the  precession 
magnitude.  This  strong  dependence  can  be  characterized  by  the  inequality 


N  »G 


whereas  as  majority  of  conventional  oscillators  is  characterized  by  the  opposite  inequality. 

It  is  clear,  that  such  strong  nonlinearity  of  STO  leads  to  a  significant  broadening  of  the 
generation  linewidth.  Therefore,  for  practical  purposes  it  would  be  desirable  to  work  in  the 
region  where  the  nonlinear  frequency  shift  N  is  as  small  as  possible. 

It  is  known  that  the  nonlinear  frequency  shift  N of  magnetization  precession  strongly  depends 
on  the  orientation  of  magnetization  of  the  magnetic  layer.  This  effect  allows  one  to  control  the 
coherent  properties  of  self-sustained  magnetization  precession  and  to  choose  optimal 
conditions  for  STO  operation. 


Nonlinear  frequency  shift  N!2x  (GHz) 


Angular  dependence  of  the  linewidth  17 
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Nonlinear  frequency  shift  coefficient  N strongly  depends  on  the 

orientation  of  the  bias  magnetic  field 
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For  an  isotropic  magnetic  film  nonlinear  frequency 
shift  coefficient  N depends  on  the  out-of-plane 
angle  O0  that  the  external  bias  magnetic  field  HO 
makes  with  the  film  plane.  In  two  limiting  cases  of 
in-plane  (O0  =  0)  and  normal  (#0  =  nil) 
magnetizations  expressions  for  N  have  the  form 

N(0  q)  2wM0J1T0  4nM0  +  4H0 
Vo  '  ^AkM0  +  H0  4 kMq  +  2 H0 

N(90  =  n/2)  = 

Here  H0  is  the  bias  magnetic  field. 

For  intermediate  magnetization  angles  N can  be 
found  only  numerically.  The  dependence  of  TVon  &0 
is  shown  in  the  figure.  It  is  clear  that  there  is  a 
magnetization  angle  0O  «  80  deg,  at  which  N 
vanishes  and  STO  does  not  experience  nonlinear 
linewidth  broadening.  At  this  point  one  expects  to 
observe  the  smallest  STO  linewidth. 


Dependence  of  the  nonlinear  frequency  shift  N 
of  isotropic  magnetic  film  on  the  out-of-plane 
magnetization  angle  O0.  Saturation 
magnetization  4zrM0  =  8  kG,  magnitude  of  the 
bias  magnetic  field  H0  =  10  kOe. 

Note  that  N changes  sign. 
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Such  angular  dependence  of  STO  linewidth  has 
been  observed  in  the  recent  experiments  [W.H. 
Rippard  et  al.,  Phys.  Rev.  B  74,  224409  (2006)] 
(see  figure). 

As  it  is  clear  from  the  figure,  the  nonlinear 
frequency  shift  coefficient  N determines  both 
the  magnitude  of  the  STO  generation  linewidth 
and  its  qualitative  angular  dependence. 

The  minimum  linewidth  is  achieved  at  the 
magnetization  angle  0O  »  80  deg,  where  the 
nonlinear  frequency  shift  N  vanishes. 

However,  almost  normal  magnetization  (0O  « 

80  deg)  requires  use  of  large  bias  magnetic 
fields  H0  ~  4 7rM0  to  saturate  the  free  layer  of 
STO.  Thus,  another  way  of  reduction  of 
influence  of  the  nonlinear  frequency  shift  N on 
the  generation  linewidth  is  desirable. 


Dependence  of  the  STO  generation  linewidth 
on  the  out-of-plane  magnetization  angle  00- 
Points  -  experiment  [W.H.  Rippard  etal.,  Phys. 
Rev.  B  74,  224409  (2006)],  red  line  - 
theoretical  dependence,  blue  line -theory  for 
a  linear  oscillator  without  nonlinear  frequency 
shift  (N=  0),  multiplied  by  10. 
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In  this  work,  we  consider  generation  linewidth  of  STO,  based  on  anisotropic  in-plane 
magnetized  "free"  magnetic  layer.  It  is  known,  that  even  small  anisotropy  strongly  modifies 
nonlinear  properties  of  magnetization  precession  and  can  lead  to  a  significant  dependence  of 
the  nonlinear  frequency  shift  coefficient  N on  the  in-plane  magnetization  angle  that  the  in¬ 
plane  bias  magnetic  field  H0  makes  with  the  easy  anisotropy  axis  of  the  layer. 

At  the  same  time,  magnitude  of  the  bias  magnetic  field  H0,  necessary  to  saturate  in-plane 
magnetized  magnetic  film,  is  substantially  smaller  than  4 7rM0  and  is  of  the  order  of  anisotropy 
magnetic  field  HA  «  4 nM^. 

An  additional  advantage  of  the  in-plane  magnetized  geometry  is  that,  for  the  same  value  of  the 
bias  magnetic  field  H0,  the  ferromagnetic  resonance  frequency  co0  is  substantially  higher  than 
for  the  out-of-plane  magnetized  case. 


In  our  approach  the  effective  anisotropy  field  HA  accounts  both  for  the  crystallographic 
anisotropy  of  the  "free"  layer  and  for  the  shape  anisotropy  caused,  in  the  case  of  STO  based  on 
magnetic  nano-pillar,  by  non-circular  shape  of  the  pillar. 


Parameters  of  anisotropic  film 
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To  find  the  generation  linewidth  of  STO,  we  need  to  obtain  expressions  for  the  damping  rate  ro 
and  nonlinearity  coefficient  N of  in-plane  magnetized  magnetic  layer.  We  derived  analytical 
expressions  for  these  coefficients: 


r0  -  ccQA 


N  = 


co0A 


(2 A2  -  3 AB  +  B2)b  +  coa  (2  A2  +  B2  )cos(2^) 


+ 


3  oj 


4  co, 


aG  is  the  Gilbert  damping  constant,  coQ  =  ^](coH  -  sin2  $  coA \coH 


(A  +  B^2A2  +2AB  +  B1  )sin2  (20) 


+  coM )  is  the  FMR  frequency, 


A  -  coH 


+  \co 


M 


+  sin 


coH  =yH  au=  AnyM  0  coA  =  yH  A 

and  //and  ^are  the  magnitude  and  in-plane  angle  of  the  internal  magnetic  field,  which  are 
connected  with  the  external  values  H0  and  by 


H  cos  (/>  =  H0  cos  4  +  H A  cos  (/) 


H  sin  ^  =  H0  sin 


Frequency  and  damping  rate 
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Analysis  shows  that  both  the  FMR  frequency 
co0  and  Gilbert  damping  rate  ro  only  weakly 
depend  on  the  in-plane  magnetization  angle 
(see  figure).  Thus,  in  the  most  interesting 
case  of  small  magnetic  field  H0  «  4 nM0 
expressions  for  the  FMR  frequency  for 
magnetizations  along  the  easy  (^0  =  0)  and 
hard  (^0  =  nil)  axes  have  the  form 

<sy  *  r^M0(H0  +  HA) 

^  *  yj4xM0(H0-HA) 

and  differ  by  only  30  %  for  H0  =  4  HA  (see 
figure). 

Gilbert  damping  rate  ro,  which  determines 
both  the  threshold  current  7th  and  range  of 
variations  of  the  generation  linewidth  lAco,  is 
almost  independent  of  </>Q,  due  to  a  strong 
ellipticity  of  magnetization  precession  in  in¬ 
plane  magnetized  film:  ro  «  aGcoJl. 


Dependence  of  the  FMR  frequency  co0  (red  line,  left 
axis)  and  Gilbert  damping  rate  r0  (blue  line,  right  axis) 
on  the  in-plane  magnetization  angle  ^0.  Bias  magnetic 
field  H0=  1.2  kOe,  saturation  magnetization  4 nM0  =  8 
kG,  anisotropy  field  HA  =  0.3  kOe,  Gilbert  damping 
parameter  aG  =  0.01. 
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The  nonlinear  frequency  shift  coefficient  TV 
strongly  depends  on  the  magnetization  angle 
(see  figure).  In  the  limit  H0  «  4 nM0  one  can 
obtain  simple  expressions  for  TV for  magnetizations 
along  the  easy  (^0  =  0)  and  hard  (^0  =  nil)  axes: 

a-r(H0-4HA)^M_ 

2cOo 

For  the  easy-axis  case  (^0  =  0)  TV  is  always  negative, 
whereas  for  =  nil  and  magnetic  field  in  the 
range  HA<H0<  4 HA  nonlinear  frequency  shift  TV  is 
positive  (see  curve  1).  In  this  field  range  TV 
vanishes  for  a  certain  intermediate  magnetization 
angle,  where  the  minimum  linewidth  of  STO 
should  be  observed.  For  larger  fields  |TV| 
monotonically  decreases  when  rotates  from  the 
easy-  to  hard-axis  orientation  (see  curves  2  and  3) 
and  STO  generation  linewidth  should  have 
minimum  for  the  hard-axis  magnetization  (^0  = 


Dependence  of  the  nonlinear  frequency  shift  N 
on  the  in-plane  magnetization  angle  for 
several  values  of  the  bias  magnetic  field  H0:  1  - 
0.6  kOe,  2-1.2  kOe,  3-1.8  kOe.  Other 
parameters  of  magnetic  film:  4 nM0  =  8  kG, 
anisotropy  field  HA  =  0.3  kOe. 
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In  this  figure  we  showed  the  calculated 
dependence  of  the  STO  generation  linewidth 
2A&>  on  the  in-plane  magnetization  angle  for 
several  values  of  the  bias  field  //0.  One  can  see 
that  this  dependence  virtually  reproduces  the 
angular  dependence  of  TV2,  since  all  other 
parameters  of  STO  depends  on  only  weakly. 

In  particular,  for  small  bias  field  H0  <  4 HA  (see 
curve  1)  2A co  has  a  strong  minimum  for  a  certain 
intermediate  magnetization  angle,  whereas  for 
large  bias  fields  H0  >  4 HA  (see  curves  2  and  3)  the 
generation  linewidth  2A&>  monotonically 
decreases  with  the  increase  of  and  has  a 
minimum  for  the  hard-axis  magnetization  = 
7t!2. 


Dependence  of  the  generation  linewidth  2A&>of 
STO  on  the  in-plane  magnetization  angle  for 
several  values  of  the  bias  magnetic  field  H0:  1  - 
0.6  kOe,  2-1.2  kOe,  3-1.8  kOe.  Other 
parameters  of  STO:  4 7rM0  =  8  kG,  HA  =  0.3  kOe, 
aG  =  0.01,  Q  =  3,  "free"  layer  thickness  d=  3  nm, 
shape  of  the  nano-pillar  -  circular  with  the  radius 
R  =  50  nm,  spin-polarization  efficiency  s=  0.2, 
bias  current /=  3  mA,  temperature  T=  300  K. 
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Recently,  experimental  study  of  generation 
linewidth  of  in-plane  magnetized  anisotropic 
STO  has  been  performed  in  [K.  V.  Thadani  et 
al.,  arXiv:0803.2871  (2008)].  In  this  figure  we 
showed  comparison  of  our  theoretical  results 
(red  line)  with  experimental  data  (dots)  for 
exchange-biased  STO  (see  Fig.  2c  in  [K.  V. 
Thadani  et  al.,  arXiv:0803.2871  (2008)]).  One 
can  see  good  qualitative  and  quantitative 
agreement  between  the  experiment  and  our 
nonlinear  theory. 

At  the  same  time,  the  linear  theory,  that  does 
not  take  into  account  nonlinear  frequency  shift 
N (see  blue  line  in  the  figure),  predicts  much 
smaller  quantitative  values  of  the  generation 
linewidth  2A&>and  much  smoother 
dependence  of  the  linewidth  on  the 
magnetization  angle  0O. 


Comparison  of  experiment  (dots)  from  Fig.  2c  in  [K. 
V.  Thadani  et  al.,  arXiv:0803.2871  (2008)]  with  our 
theory  (red  line).  Parameters:  4 7rM0  =  8  kG,  H0  = 
1.08  kOe,  Ha  =  0.2  kOe,  crG  =  0.015,  Q  =  3,  "free" 
layer  thickness  d  =  4  nm,  shape  of  the  nano-pillar  - 
elliptical  150  nm  x  50  nm,  spin-polarization 
efficiency  e=  0.32,  bias  current /=  5  mA, 
temperature  T  =  300  K.  Blue  line  shows  theory  for 
a  linear  oscillator  [N=  0),  multiplied  by  10. 
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In  this  figure  we  showed  comparison  of  our 
theory  (red  line)  with  experimental  data  (dots) 
for  thick-fixed-layer  STO  (see  Fig.  2d  in  [K.  V. 
Thadani  et  al.,  arXiv:0803.2871  (2008)]).  Again, 
the  agreement  between  the  experiment  and 
our  nonlinear  theory  is  rather  good. 

In  both  experimental  cases  the  bias  magnetic 
field  H0  was  larger  than  then  4 HA  and, 
according  to  our  theory,  generation  linewidth 
of  STO  2A co  monotonically  decreases  with  the 
increase  of  the  magnetization  angle  and  has 
a  strong  minimum  for  the  hard-axis 
magnetization  =  n!2. 


Comparison  of  experiment  from  Fig.  2d  in  [K.  V. 
Thadani  etal.,  arXiv:0803.2871  (2008)]  with  our 
theory  (red  line).  Parameters:  4 7rM0  =  8  kG,  H0  = 
1.2  kOe,  HA  =  Q.l  kOe,  crG  =  0.015,  Q  =  3,  "free" 
layer  thickness  d  =  4  nm,  shape  of  the  nano-pillar 
-  elliptical  130  nm  x  70  nm,  spin-polarization 
efficiency  s=  0.375,  bias  current  1=  5  mA, 
temperature  T  =  300  K.  Blue  line  shows  theory 
for  a  linear  oscillator  (N=  0),  multiplied  by  10. 
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Our  theoretical  results  suggest  an  optimum  design  of  STO,  in  which  the  influence  of  thermal 
fluctuations  will  be  minimized  and  the  generation  linewidth  2A&>will  have  minimum  possible 
value. 

To  achieve  this  goal,  it  is  necessary  to  operate  STO  in  a  regime,  when  the  nonlinear  frequency 
shift  TV  vanishes  and  the  generation  linewidth  does  not  experience  nonlinear  broadening. 

The  simplest  way  to  do  this  is  to  use  in-plane  magnetized  along  the  hard  axis  nano-pillar  with 
the  bias  magnetic  field 


H0  =  4  Ha 


The  generation  frequency  of  STO  in  this  case  is  equal  approximately  to 


and  can  be  tuned  by  changing  the  anisotropy  field  HA  (e.g.,  by  changing  the  aspect  ratio  of 
ellipsoidal  magnetic  nano-pillar  forming  "free"  layer  of  STO). 
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•  We  have  developed  theory  of  generation  linewidth  of  spin-torque  oscillator  (STO) 
with  an  in-plane  magnetized  anisotropic  "free"  magnetic  layer. 


•  The  dependence  of  the  generation  linewidth  of  STO  on  the  in-plane  magnetization 
angle  is  determined  mostly  by  the  angular  dependence  of  the  nonlinear  frequency 
shift  coefficient. 


•  For  relatively  small  bias  magnetic  fields  (HA  <  H, o  <  4 Ha)  the  generation  linewidth  has 
a  minimum  for  a  certain  intermediate  magnetization  angle,  whereas  for  large  bias 
fields  (H0  >  4 Ha)  the  linewidth  has  a  minimum  value  for  the  magnetization  direction 
along  the  hard  in-plane  axis. 

•  Our  theory  provides  qualitatively  and  quantitatively  correct  description  of  recent 
experimental  studies  of  the  angular  dependence  of  the  generation  linewidth  in  in¬ 
plane  magnetized  anisotropic  STO. 

•  Our  results  allows  one  to  select  an  optimum  design  of  STO,  in  which  generation 
linewidth  will  have  a  minimum  possible  value. 


